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Abstract. The paper deals with the asymptotic behavior of solutions to a non-local diffusion 
equation, ut = J *u — u ~ Lu, in an exterior domain, Q, which excludes one or several holes, 
and with zero Dirichlet data on \n. When the space dimension is three or more this behavior 
is given by a multiple of the fundamental solution of the heat equation away from the holes. 
On the other hand, if the solution is scaled according to its decay factor, close to the holes it 
behaves like a function that is L-harmonic, Lu = 0, in the exterior domain and vanishes in 
its complement. The height of such a function at infinity is determined through a matching 
procedure with the multiple of the fundamental solution of the heat equation representing the 
outer behavior. The inner and the outer behavior can be presented in a unified way through a 
suitable global approximation. 



1. Introduction 

Let ^ C M^, > 3, be a bounded open set with smooth boundary and let Q = \ H. 
We do not assume H = i}'' to be connected, so that it may represent one or several holes in an 
otherwise homogeneous medium. Our goal is to study the large-time behavior of the solution to 
the nonlocal heat equation in that exterior domain with zero data on the boundary, namely. 



(LI) 



dtu{x,t) = / J{x — y)(u{y,t) — u{x,t)) dy in r2x(0, oo), 
Jrn 

u{x,t) = in ^ X (0, oo), 

u{x,0)=uo{x) in M^, 

with a kernel J that is assumed to be a nonnegative continuous function with unit integral. 

We will restrict ourselves to the case where J is smooth, radially symmetric, with a compact 
support contained in the unit ball centered at the origin and J(0) > 0. However, results similar 
to the ones we obtain here should hold for more general kernels. 

The hypotheses on the initial data, mq, are that they are nonnegative, bounded, integrable, 
and identically zero in the hole 7i. We also assume, without loss of generality, that £ T-L. 

Evolution problems with this type of diffusion have been widely considered in the literature, 
since they can be used to model the dispersal of a species by taking into account long-range 
effects, [1], [7], [H]. It has been also proposed as a model for phase transitions, [2], [3], and, 
quite recently, for image enhancement, [lOj . 
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The nonlocal equation in (jl.ip has been the subject of many recent works in the case where 
the spatial domain is M.^ , and also when it is a smooth bounded domain and it is complemented 
by a suitable Dirichlet or Neumann type boundary condition, see the recent monograph [1] and 
the references therein. However, up to our knowledge, this is the first time that an exterior 
domain is being considered. 

When there are no holes, ?^ = 0, mass is conserved, J^m u{x,t) dx = f^MUo{x)dx, and 

the solution to behaves for large times as the solution v to the local heat equation with 
diffusivity 

(1.2) «^=^/ \z\^J{z)dz 

and the same initial data [8], [13]. More precisely, 

(1.3) lim 

max \u{x, t) — v{x, t)\ =0. 



t—>oo 



Hence the asymptotic behavior of u can be described in terms of the fundamental solution of 
the heat equation with diffusivity a, Fq. This solution has a self-similar structure. 



TUx,t) = t-'^I^U^ (^) , U^{y) = (47ra)-^/2 



e 



Therefore, in self-similar variables we have convergence towards the stationary profile MUa, 
where M = J^m uq{x) dx, that is, 

(1.4) lim max |t^/2^(yt^/^ t) - MC/«(y)| = 0. 

Thus, there is an asymptotic symmetrization: no matter wether the initial datum is radial or 
not, the large time behavior of u is given by a radial profile, which, of course, has the same mass 
as the datum. 

In the presence of holes the situation is very different. On one hand, mass is not conserved. 
On the other hand, the presence of the hole breaks (in general) the symmetry of the spatial 
domain, and an asymptotic symmetrization is no longer possible. Nevertheless, when > 3 we 
still have a conservation law, J^^ u{x, t)(j){x) dx = J^n uo{x)(j){x) dx, with cp the unique solution 
to 

J ^ (j) = (f) in 17, 

(1.5) \(t> = Q in n, 



>(x)^l 



as X — )• oo. 



Moreover, there is a non-trivial asymptotic mass which coincides with the conserved quantity 
M* = J^ff uo{x)<j){x) dx. Besides, if we stand far away from the holes, they will be seen as a 
point, and we may still expect some symmetrization. Indeed, we will prove that 

lim t^/'^u{yt^/'^,t) = M*Ua{y) uniformly in \y\>5> 0. 

t— >oo 

The only effect of the holes in this outer limit is the loss of mass. However, close to the holes, 
in the inner limit, their effect is much more important. On compact sets solutions still decay as 
Q^^-N/2y jf -^g scale the solutions accordingly, we get that the new variable w{x, t) = t^^'^u{x, t) 
satisfies 

OtW = J * w — w ^ . 

2 t 
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Thus, we expect w to converge to a L-harmonic function $ with zero 'boundary data', 

L$ := J* ^) - $ = 0, x^n, $ = 0, x^n. 

This problem does not have uniqueness. However, a unique solution can be determined by 
prescribing the value at infinity. Indeed, if $ — )• C* as — )■ oo, then ^ = C*(j). 

In order to select the right constant C* characterizing the asymptotic limit of u;, we need some 
extra information. This will come from the outer limit, following a typical matched asymptotics 
procedure. Let us explain in some detail how this is done. 

In view of the asymptotic behavior of 0, we can describe the outer limit in the alternative 
form 

t^/'^\u{x,t) - M*(f>{x)Ta{x,t)\ ^ as t oo uniformly in > (^t^/^, 5 > 0. 

Besides, since lim t^/'^Ta{x, t) = (47ra)~^/^ uniformly on compact sets, the expected inner limit 
can be written as 

as t — 7- oo uniformly in compact sets. 

We will prove that we can go beyond compact sets, and that this limit holds in a set of the form 
\x\ < for some 6 > 0. Hence, there is an overlapping region between the inner and the 

outer development, and they can be matched, leading to 

C* = M*(47ra)-^/2_ 

Notice that we have been able to describe the inner and the outer behavior in a unified way. 
Hence we may gather the two results in just one theorem. 

Theorem 1.1. Let N > 3, < uq e L'^{R^) n L°°(M^) and M* = J^j^ uo{x)<j){x) dx. Then, 
(1.6) t^/'^\u{x,t) - M*(l3{x)Ta{x,t)\^Q as t ^ oo uniformly in R'^ . 

Remark 1.1. The set Q + i?(0;l/2) has only one unbounded component U. This component 
l/( may contain several components of 0,. Though only one of them is really unbounded, due 
to the non-local character of the diffusion operator, all of them are 'connected to infinity'. 
However, there may exist components of which do not intersect U, and which are, thus, 
'truly' bounded. The function (j) is identically there. Therefore, the scaling we are using is 
not adequate to characterize the asymptotic behavior in such subdomains. Indeed, the decay 
rate in these kind of bounded 'components' is exponential, and the asymptotic profile is an 
eigenfunction of the operator L with zero Dirichlet boundary conditions in the complement of 
the 'component', associated to the first eigenvalue, [8]. 

Remark 1.2. The asymptotic behavior for the case where there are no holes ()1.4p can also be 
described in the form ()1.6p . since the solution to (|1.5p when Ti = 9 is (j) = 1. 

When N < 2, mass is expected to decay to zero (logarithmically when N = 2, and like a 
power when = 1), and the solutions will decay faster than 0{t~^/'^). The analysis becomes 
much more involved, and is postponed to a future paper. 

An analogous study in an exterior domain has been performed for the corresponding local 
problem in jlTj . However, the author makes extensive use in the proofs of the explicit form 
of the fundamental solution of the local heat equation in the whole space. One of the main 
difficulties in the non-local case is that we do not have such an explicit expression. 
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An alternative approach, in which fundamental solutions do not play such a key role in the 
proofs, has been used for nonlinear (local problems), see for instance [6], [12]. However, the 
operators under consideration in these cases are invariant under some scaling transformations, 
which is not the case in the present problem. 

Nevertheless, although it will not be apparent from our proofs, scaling still plays a key role 
in the outer region. In fact, as observed in [9], under the usual parabolic scaling our operator 
'converges' to aA. This is the reason why the outer behavior is given by the fundamental 
solution to the heat equation with diffusivity a. This idea has been used in [TB] to study the 
asymptotic behavior for the nonlocal heat equation with absorption in a certain critical case. 

We complete the study of the asymptotic behavior by analyzing the function at infinity. 
We prove that all the derivatives of (p behave at infinity as those of the fundamental solution of 
the Laplace operator. This result might be of independent interest. 

Notations. In what follows we will denote Lu{x,t) := J^m J{x — y)[u{y,t) — u{x,t)^ dy. 

2. Preliminaries 

On the notion of solution. Let uq G L-'^(M^), uq = a.e. in %. A solution of (jl.ip is a 
function u £ C([0,oo) : ^^(R^)) such that for alH > 

(2.1) u{-,t) = U(j+ / (^J*u{-,s) — u{-,s))ds a.e. in f], u{-,t)=0 a.e. in H. 

Jo 

Subsolutions and supersolutions are defined as usual by replacing the equalities in the definition 
above respectively by < or >. 

Remark 2.1. If n is a solution, then u G L^(M^ x [0,r]) for all T > 0. Hence, ([Ll]) holds, 
not only in the sense of distributions, but also a.e. in x (0,oo). Moreover, we also have 
u € C^([0, oo) : L^(M^)), and the equation holds a.e. in for all t > 0. 



Existence and uniqueness. Existence and uniqueness of a solution to (|l.ip can be proved 
for any uq G L^(]R^) such that uq = a.e. in Ti by using a fixed point argument. We do not 
need to assume neither uq G L°°(M^), nor a sign restriction. 

Theorem 2.1. For any uq G L^(M^) such that uq = a.e. in Ti. there exists a unique solution 

of dni). 

Proof. The proof is quite similar to the one for the case Ti = 9 performed in [5]. 

Let Bto = {n G C([0,to] : L^(M^)) : u{-,t) = a.e. in n for ah t G [0,to]}- Endowed with the 
norm 

llklll = max ||n(-,t)||ii(iRiV), 

it is a Banach space. We define the operator T : Bt^ — Bt^ through 

ft 



(2.2) {Tu){.,t) 



^o(") + / {J * u{- , s) — u{- , s)) ds a.e. in ri, 
Jo 

a.e. in H. 
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This operator turns out to be contractive if to is small enough. Indeed, 
/ \Tip-TiP\{;t) 

{\\J\\LHm) + l) ll(<^ - '^)(-, s)llii(iRiV) ds. 

Hence, IHT^^ — TV'IH < 2to|||¥' ~ V'lll- Thus, T is a contraction if to < 1/2. Existence and 
uniqueness in the time interval [0, to] now follow easily, using Banach's fixed point theorem. 
Since the length to of the existence and uniqueness time interval does not depend on the initial 
data, we may iterate the argument to extend the result to all positive times. □ 




Comparison. Comparison is an immediate consequence of the following T-contraction prop- 
erty in L^. Again, we are only assuming that the initial data are integrable. 

Theorem 2.2. Let ui and U2 be two solutions of (jl.ip with initial data ui^0iU2fi S L^(R^). 
Then, for every t > 0, 

{Ui - U2) + {-,t) < / {uifi-U2,o) + - 



Proof. We subtract the equations for ui and U2 and multiply by Ij^^^u^}- Since ui — U2 € 
CH[0,oo);Li(M^)), 

dt{ui - U2)l{u^>u2} = dt{ui - U2)+. 
On the other hand, just using that < l|„j>„2} — 1 ^^"^ J > we get, 

J *{ui- 'U2)l{„i>„2} < J*{ui- U2) + . 

Finally, {ui — U2)'^{ui>u2} = (""i ~ ^2)+- We end up with 

{J * (til — ^2)+ — {ui — U2)+ a.e. in Q, 

a.e. in T-L. 
Integrating in space, and using Fubini's Theorem, we get 

dt [ iui-U2)+i-,t) < 0. 



□ 

Remark 2.2. The same proof applies when there are no holes, Q = M^. 

As a corollary we have, on one hand comparison and, on the other hand an L^-contraction 
property for solutions, 

||(ni - n2)(-,t)||ii(ffiiV) < \\UI,0 -U2,o\\L^RNy 

The latter property implies the continuous dependence of solutions on the initial data, and a 
uniform control of the L^-norm along time, 

- II'"oIIli( 
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The proof of Theorem 12.21 apphes under much weaker hypotheses, leading to a more general 
comparison principle that will be useful in the sequel. 

Theorem 2.3. Let {Q{t)}t>o be a continuous (with respect to the HausdorjJ distance) family of 
open smooth domains o/M^. Let ui,U2 G C^{[0,oo); L^{M.^)) be such that 

dfUi — Lui < 0, dtU2 — Lu2 > for a.e. x £ Vt > 0, 

and ui{x,t) < U2{x,t) for a.e. x € \ 0,{t), \/t > 0. Then, for every t >0, 

(ui - ii2)+(-,t) < / (ni(x,0) - n2(x,0))+. 



Remark 2.3. In order to compare a subsolution and a supersolution in fl{t), t > 0, we only need 
them to be ordered in {n{t) + B{0; 1)) \ n{t), t > 0, and in n{0) + B{0; 1). 

Time decay. Solutions to problem (jl.ip satisfy 

e-^uo + e-(*-") J * n(-, s) ds a.e. in \ U, 



(2.3) u{-,t) 



a.e. in %. 



Therefore, u has the same spatial regularity as the initial datum uq, but not more. In particular, 
if no is not bounded, neither is u{-,t) for any later time. 

On the contrary, if the initial datum is bounded, the solution stays bounded as time goes by. 
Moreover, if the datum belongs to L^(]R^) n L°°(]R^), the solution decays like 0(t"^/2). To 
check this, we first observe that, since u is a solution to (jl.ip and u{-,t) = a.e. in T-L for all 
t G [0, oo), then 

(2.4) ut = Lu - X-uiJ * u) in X (0,oo). 

Besides, if uq > 0, we have u{x,t) > 0. Hence, u, which belongs to C([0, oo); L^(M^)), is a 
subsolution to the Cauchy problem (no holes). Therefore, it lies below the solution uc to the 
Cauchy problem with the same initial datum. 

On the other hand, as we have already mentioned (see dL!])), if tio G L^(M^) n L°°(M^), then 
uc = 0(r^/2). Thus, 

(2.5) 0<u{x,t) <Ct-^/'^. 

Remark 2.4. In order to obtain the asymptotic behavior of uq, the paper [8] requires uo,uo S 
L^(IR^). The paper [13] dispenses with the integrability assumption on uq, and on the sole 
hypothesis uq G L^{R^) proves that 

lim t^/^ max \u{x,t) — e~^uo{x) — v{x,t)\ = 0, 



where v has the same meaning as in (|1.3p . Hence, the hypotheses uq G L^(M^) n L°^(M^) are 
enough to have (|1.4p . 



L°°-SOLUTiONS. Instead of the L^-theory outhned above, an L°°-theory for (jl.ip is also 
possible. Given an initial datum uq G L°°(]R^), an L°°-solution to problem (jl.ip is a function 
u G C([0,oo);L°°(]R^)) such that ([23]) holds for ah t > 0. 
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As in the case of L^-solutions, existence and uniqueness follow from a fixed-point argument. 
Indeed, Bto = {u £ L°°(M^ x [0,to]) : u{-,t) = a.e. in n for all t G [0,to]} endowed with the 
norm 

llklll = ^max ||u(-,t)||i«=(]K]V), 

is a Banach space, and the operator T : Bt^ — )■ Bt^ defined by ()2.2p satisfies \\\Tip — TtP\\\ < 
2to|||¥' ~ V'lll- Hence, it is contractive if to < 1/2, which implies the local existence and unique- 
ness result. By iteration, taking as initial datum u{-,to) G L°°(R^), we obtain existence and 
uniqueness for [0, 2to] and therefore for all times. 

Besides, L°°-solutions depend continuously on the initial data. Indeed, let ui and U2 be L°°- 
solutions with initial data ui^q, U2fl G L°°(]R^). Since Uj satisfies (|2.3p (with uq = Ui^), an easy 
computation shows that for every t > 0, 

max \\ui{-,s) - U2{-,s)\\Lac(j^N) < \\ui^o - U2,q\\l°°(rn). 

In particular, the maximum principle for bounded solutions applies, namely 

\H-:'t)\\L^{R'^) < lh(->0)||ioo(iR>iV) \ft > 0. 

Essentially the same computation yields a comparison result for bounded sub- and superso- 
lutions if the hole does not change with time. Let us see this in some detail. Assume ui^ < U2fi 
in R^. Let ui be a bounded subsolution with initial datum ui^ and U2 a bounded supersolution 
with initial datum U2fi in Q, = \ H with ui < U2 in "H for every t > 0. Let u = ui — U2 and 

= ^^1,0 - ""2,0- Then, 



e~*no(x) / e"^*~*^ J * ^(x, s) if x G r2, 

Jo 

iixen. 



u{x,t) < < 



Hence, for x G O, 



with leads to 



t 

u{x,t) < I e~^^~''^J*u{x,s)ds< supn(l-e"*), 

M^x(0,t) 

sup u < sup u {1 — e~*). 

RJVx(0,i) RJVx(0,i) 

Therefore, supj^jv x (o,t) u < for every t > 0. This means that u < 0. Thus, ui < U2- 

Remark 2.5. The comparison result for bounded sub- and supersolutions of the Cauchy problem 
was proved in [15| using a different technique. The result of that paper is in fact more general, 
since it also applies to semilinear equations, as long as the nonlinearity is locally Lipschitz 
continuous. In addition, its proof copes without further ado with the case of holes that change 
with time. 

Remark 2.6. From now on we will always assume that < uq G L-'^(IR^) n L°^(M^). 



Representation formula. By (j2.4p . u can be expressed in terms of the fundamental solution 
F = F{x, t) to the operator dt — L in the whole space by means of the variations of constants 



8 CORTAZAR, ELGUETA, QUIROS AND WOLANSKI 

formula. Thus, for t > to we have 

{2.6) u{x,t) = F{x-y,t-to)u{y,to)dy- / F{x - y,t - s)X-H{y){J *u{-, s)){y) dy ds. 

JR^ JtQ JRJV 

The fundamental solution can be decomposed as 

(2.7) F{x, t) = e~*(5(x) + uj{x, t), 

where d{x) is the Dirac mass at the origin in and u is defined via its Fourier transform as 

(2.8) 6j{^,t) = e-'{e'^^^'-l), 

see [5], from where it is easy to see that a; is a smooth function and 

uj{x, t) dx = cj(0, t) = 1 — e~*. 

Combining (j2.6p and (j2.7p . we obtain 

n(x,t) = e"(*"*'^)n(x,to) + / (^{x - y,t - to)u{y,to) dy 

JrJv 

rt 



(2.9) - / e-'^'-'^Xn(.x){J*u{;s)){x)ds 

J to 



to 

/ ^{x - y,t - s)X-uiy)iJ * u{-,s)){y)dyds. 

to Jm.^ 

This formula will be the starting point to obtain both the behavior at infinity of <j) and the outer 
large time behavior of u. 

Estimates for w. In order to take profit of the representation formula ()2.9p . we need good 
estimates for the regular part, w, of the fundamental solution. 

A first source for estimates comes from the asymptotic convergence of uj to the fundamental 
solution of the local heat equation with diffusivity a, which is proved using Fourier transform 
techniques [13]. Indeed, given a multi-index /3 = (/3i, . . . ,/3Ar) G Nq of order = s, 

(2.10) t(^+')/2||L>'^a;(-,t)-L»^r„(-,t)||ioo(iRiV) ^0 as t ^ oo, 
where D^f = dxl ■ ■ ■ dx^ f. In particular, 

(2.11) \D(^uj{x,t)\ < Ct^^ if |/3| = s. 

These estimates give the right order of time decay. However, they do not take into account the 
spatial structure of oj, and will not be enough for our purposes. We need to know something 
about the spatial decay of a; as \x\ — t- oo. 

A second source of estimates is the expansion 

rj*"(x) 



oo 



(2.12) uj{x,t) = e-^Y^ 



n=l 



which follows by using the Taylor series of the exponential in (j2.8p . This expression, recently 
derived in [5], was used by the authors to obtain estimates giving the behavior of a; as oo 
for each fixed t. However, though they give the right order of spatial decay for each time, they 
become very poor as t — )■ oo. 
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Hence we need a different approach: we will obtain 'good' estimates through a comparison 
argument. Indeed, as observed in [16], is a solution to 

,^ \dtu:{x, t) - Luj{x, t) = e-V(x) in x (0, 00), 

^ ' ' \w(x,0) = in R^, 

a fact that can be checked, for example, either by differentiating uj{£,, t) with respect to t or 
using that e~^6{x) +u}{x, t) is a solution of dtF — LF = 0. Good estimates will then follow from 
the use of appropriate barriers. 

Proposition 2.1. The function uj is nonnegative. Moreover, given a multi-index /3 G Nq'^ of 
order |/?| = s, 

(2.14) \Df^u{x,t)\<Cj-J^, 

(2.15) / \D^uj{x,t)\dx <Ct-'/^. 

The cases s = 0, 1 were already proved in |16j . through a clever use of the above mentioned 
comparison principle from [15]. The proof of the general case is completely analogous. The fact 
that UJ is nonnegative was also proved in [16] and can also be seen from the series expansion (j2.12p . 

Remark 2.7. Estimates (|2.14|) and (|2.15|) are in some sense optimal, since they are invariant 
under the 'parabolic' scaling {Tf^u){x,t) = k^u{kx,k'^t). This scaling will play a major part in 
some of our proofs. 



3. The stationary problem 



As we have seen in the introduction, the L-harmonic function (j) solving ()1.5p plays a key 
role in the description of the asymptotic behavior of solutions to (jl.ip . In this section we prove 
that problem (|1.5p is well posed, and we obtain some properties of (p that will be required 
later. Though not needed in the sequel, we complete our study of the stationary solution (p by 
describing precisely its behavior at infinity. 



Remark on the notion of solution. In principle we only need (j) G ^^^^.(IR ), so that the 
convolution makes sense. But the equation implies that (j) has in Q the same regularity as J. 
Hence <j) is at least continuous in and, in our case, (j) is smooth in Q. However, (p is not 
continuous across the boundary of fi. 



3.1. Existence and uniqueness. To prove existence we will approximate the domain Q by a 
sequence of bounded domains, = H i3(0; n). Let us start by studying the Dirichlet problem 
for such bounded domains. 

Lemma 3.1. Let G C R^ open and bounded. 

(i) Given f € L°°(G^), there exists a solution to problem 

L(p{x) = in G, (j) = f{x) in G'^. 
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(ii) Let 01, (/>2 G L°°(M^) n C(G) such that 

Lcpi > L(j)2 in G, (pi < (j)2 in . 

Then 0i < (1)2 in M^. 

Proof. Existence. Let K, = {ip £ C{G) : \\4'\\l°°{g) ^ II/IIl°°(G<=)}- We define an operator 
T : /C — ?■ /C through the formula 



{Ti;){x)= / J{x - y)ij{y) dy + / J{x - y)f{y) dy. 
Jo 

It is easy to see that ||T'V'||ioc(G) < II/IIl°°(G'^) and ||T'V'||c71(g) ^ C'jII/IIl°°(G'=)- Hence, by 
Schauder's fixed point theorem, T has a fixed point (j) £ fC. Then, the function 

(i>{x), X e G, 

J{x), xeG^. 
is the sought solution. 
Comparison. Let (p = (pi — (p2- Then, 

(p{x) < / J{x-y)(p{y)dy. 
JG 

Assume for contradiction that M = max^c/) > 0. Let xq £ G he a point at which the maximum 
is attained and let Gi be the corresponding connected component. The set K = {x £ Gi : 
<p{x) = M} is nonempty and closed. On the other hand. 



M = ^{xq) < / J{xo - y)(P{y) dy < M J{xq -y)dy< M. 
JG JG 

This is only possible if (p{x) = M in G n supp J{xq — •). Since J(0) > 0, this implies that 
(p{x) = M in G n B{x{)]5) for some (5 > 0; hence K is an open subset of Gi. Since Gi is 
connected, we deduce that K = Gi., so that, (p = M > Q m. Gi. But this implies that 



M = (/)(xo) < M / J(xo -y)dy < M 
JG 

for any xq G dGi, a contradiction at any point where the density of is positive. □ 

In order to prove that the function (p obtained in the limit procedure approaches 1 at infinity, 
we will compare the approximate solutions defined in bounded sets with a suitable barrier z. 
Hence we need to estimate Lz. 

Lemma 3.2. Let z{x) = + b)~'^ , 7 > 0. There exists a number b = b{J,^,N) > such 
that 

(3.1) Lz{x) < -2a7(|x|2 + 6)-(T+i)(Af - 2 - 27), 

where a is the constant defined in (jl.2p . 



Proof. After some cumbersome, but easy, computations, using Taylor's expansion and the radial 
symmetry of J, we find that 

Lz{x) = - 2aj{\x\^ + by^^+^'>{N - 2 - 27) - 407(7 + l)(kP + b^^^+^h 
+ 0((|x|2 + 6)-(T'+2)), 
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from where the result fohows just choosing b large enough. □ 

Proposition 3.1. There exists a unique solution to (|1.5|) . Moreover, < (p < 1 and there exists 
a positive constant K such that 

(3.2) i-m< ^(^-2)72 - 

Proof. Existence. Let no in N such that H C i?(0;no). Consider the sequence of functions 
{4>n{x)}n>no-, where is the solution to 



(3.3) 



L0„(x) = O inS(0;n)\^, 

0n = in H, 

(/)„ = ! in 5(0; n)^ 



A simple comparison argument proves that < ipnix) < 1 for all n > no, x € M^. Hence, 
<t>n{x) > (j)n+i{x) in {B{0;n)y, which yields, again by comparison, < (j)n+i{x) < (j)n{x) < 1 for 
all n > no, X € R^. Therefore, the monotone limit 

(pix) = lim (pnix) 

n— >oo 

is well defined, and satisfies < (f) < 1. Letting n — )■ oo in (j3.3p . we get 

L(p{x) = in n, (p = Om U. 



We still have to prove ()3.2p . which implies in particular that 4i{x) — )■ 1 as — )■ oo. This is 
done by comparison with the barrier z described in Lemma 13.21 with 27 = N — 2. Indeed, by 
taking K large enough we have, 

i^z > 1 > in -H, i^z >!-(/.„ = in (5(0; n))^ L{Kz) < = L(l-0„) in B{0;n)\'H. 

Hence, Kz > 1 — in and, passing to the limit, we obtain ()3.2p . 

Uniqueness. Let (j) and (p be two solutions to (jl.Sp . We consider the family of functions 
(pxix) '■= <p{x) + A. Since both (p and (j) fire bounded, there is a value Ai such that (j)^-^ > in 
M^. Let 

Ao = inf{A : > in M^}. 

We shall prove that Aq = 0, which means that cp > (p in M^, from where uniqueness follows, 
since the roles of (j) and (p can be interchanged. 

Assume for contradiction that Ao > 0. Since (p{x) and (p{x) tend to 1 as |x| — )■ 00, there exists 
Ri such that (p\f^{x) > (f){x) if > Ri. Moreover, cp^^^^ > (p inH. Hence, as both (p and (f) are 
continuous in R^ \ %, there exists xq in 5(0; Ri)\'H such that 4>Xg{xo) = (p{xo). The set 

S = {xe 5(0; i?i) \ n : ^;,,(x) = (P{x)} 

is closed and non-empty. If xi is a point in dS where the density of S'^ is positive, then 
= 4'{xi) but J*(?^);^^(xl) > J*(/>(xi). Since this is a contradiction we get that Ao = 0. □ 
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3.2. Asymptotic properties. We now obtain the decay at infinity of the derivatives of (f). In 
the sequel we let = 1 — cj). 

Proposition 3.2. Given s G No, there exists a constant Cg such that, for every multi-index 
I3£N^ of order \/3\ = s, 

(3.4) \dP^(^x)\<j-^^^, xen. 

Proof. The function (j) is the unique bounded solution to the evolution problem (jl.ip with initial 
data <j). Hence, using the representation formula (|2.9p with tQ = and n(x,to) = 4'{^)j we get, 

4>{x) = . ^ -t I - y,t)(l){y)dy - X'u{x){J * 4>){x) 



1 -e- 
1 



/ uj{x-y,t- s){J*(j)){y)dyds. 
Jn 



Therefore, since J^n u^{x,t) dx = 1 — e *, we have 

i^ix) = z — - — 7 / uj{x - y,t)'il){y)dy + X'h{x){J * (t)){x) 

(3-5) ^ 

+ - -/ / uj{x-y,t- s){J*<j)){y)dyds. 

1 - e ^ Jo Jn 

Given k > 0, we denote 

^'^(x) = k^-^4>{kx), i^'^ix) = k^-^i;{kx), Jk{x) = k^J{kx), ujk{x, t) = k^{kx, kH), 
and Tik = G : kx G H}. If we scale equation ()3.5p according to these recipes, we get 

V''(x) = :-^ / u;kix-y,k-H)iP''iy)dy + Xndx){Jk*cP'')ix) 



+ , ^ „, / / ujk{x-y,k ^{t- s)){Jk*(j)'')iy)dyds. 
lo J-Hk 



l-e 



Let 5 > be given. There exists a value ks such that X-}if,{x) = if |x| > 6 and k > k^. 
Therefore, 



dP^\x) = / D^ukix - y, k-H)^l:\y) dy 



A 



1 



+ -^_^-t I I D^Mx-y,k~\t-s))iJk*(l)''){y)dyds, 



JHk 



if \x\ > 6 and k > ks- 

We start by estimating B. If y G Tik and > 5, then |x — y| > 6/2 if k > ks- On the other 
hand, since < (/>''(x) < k^~^, we also have < Jk*())'' < k^'^. So, recalling that = Ck'^ 
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and the bounds in (j2.14p - that are scahng invariant - we get 

< C5,N-, — / {t-s)ds = Cs,N-, —t'^ < Cs,N if t = 

i — e Jq i — e 

Now we estimate A with t = k"^. Recall that < ip{x) < C/|x|^~^. Therefore we have, using 
again (l2H]l and (|2J5|) . 

1^1 < [ \D^ujk{x - y, l)\r-^dy + ^''^ / - y, 1)| dy 

l-e-" J\y\<S/2 \y\ l-e-" J\y\>5/2 

Cs,N f dy Cs,N ^ 

- h\<m 

In conclusion, given 6 > there exists ks such that, if k > ks and \x\ > 6, then < 
C<5,Af- 

Let us take = 1, Ci = Ci^n, y = kx so that fc = \y\. Then, we get 

\yf-'^+'\D>^iP{y)\ = k^~'^+'\D^iP{kx)\ = \D^iP''{x)\ <Ci if \y\ = k > ki. 

This proves estimate (|3.4p except for a bounded set. However, cp is smooth in fi, and hence the 
estimate is obviously true for bounded sets. □ 

Remark 3.1. The proof of the bound ()3.4p only requires J G C"*"*"^. 

Though the second derivatives of ip decay like 0(|x|~^) as \x\ — )• oo, a special combination, 
the Laplacian, decays faster. In particular, Ai/j is integrable at infinity. 

Corollary 3.1. There exists a constant C such that \/\il){x)\ < , ,^4.2 in Q.. 

Proof. Using Taylor's expansion we get L'ip{x) = aA'0(3;) + 0( max|^l=4 ||Z)^'i/'||2,oo(5j(3,))) . Since 
ip is L-harmonic, the result follows from the estimate ()3.4p for the fourth order derivatives. □ 

Remark 3.2. The proof of Corollarv 13.11 only requires J £ C^. If we weaken this assumption 
and only require J G C^, we get |AV'(x)| < C /\x\^~^^ , and the Laplacian of il) is still integrable 
at infinity. 

We can take profit of the integrability of Ai/^ at infinity to obtain the precise asymptotic 
behavior of ijj. 

Proposition 3.3. There exists a constant C* such that for every multi-index /5 G of order 
\fi\=s, 

(3.6) \x\^-'^+' D^ij{x)-Dl^^*{x) as |x| ^ 00, i^* {x) = C*\x\^-^ . 

Proof. Scaling and compactness. From ()3.4p we get that the scaled functions -0^ satisfy 
\D^^l)^{x)\ < C/|x|^-2+'* if = s. Hence, thanks to Arzela-Ascoh's theor em, for any se- 
quence {/cnjneN there exist a subsequence, that we name again {/c^jnGNi and a function ip* such 
that, for any multi-index /?, 

(3.7) D^ip^'' D^tlj* uniformly on compact sets of \ {0}. 
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Moreover, i/j* inherits from the functions ip^ the estimate 



(3.8) 



< il^^x) < C/\x 



N-2 



Identification of the limit along subsequences. We consider the scaled operator Lk defined by 
Lf^v = k'^[Jk * V — v^. Performing a Taylor expansion, it is easy to check that Lk^p — )• aAcp 
uniformly for any G C^(M^). 

Let now ip £ C^{R^ \ {0}). We write 
(3.9) / ^P*aAip= [ il)^Lk^+ f {iP* - i;'')aAip + [ ^''{aA^ - Lk^). 

JRJV JlRjV J^N J^N 

Since L^'ip^ = 0, we have /j^jv ip^Lk^ = f^^ Lkil^^^p = 0. Besides, using ()3.4p . for k large enough 
we have 



^P^iaAif- Lkp) 



< C 



\aAip — Lf^ifl — )• as k ^ oo. 

/R^n(supp <p+Bi/fc{0)) 

Therefore, taking k = kn in (j3.9p and letting n — )• oo, we get J^^ ip*aA(p = 0, which implies 
that A^* = in \ {0}. Since |x| V* is bounded, see estimate (|3.8p . we conclude that 
there exists a constant C* such that 'il)*[x) = C* /\x\^~'^ . 

Uniqueness of the limit. We next prove that the constant C* is independent of the sequence 

{knjneN- 

A direct computation shows that 

dip* 



1 



dB{0;R) 



dr 



dS = C*{2- N) n^-^dBiO; 1)) for all R>0. 



On the other hand, if 7^ CC i?/jQ, 

f ^dS = [ ^dS+ [ AV'(x) dx for all R > Rq. 

JdB{0;R) or JdB{0;Ro) * JBrXBr^ 

Since f^^^^^ A^(x) dx has a finite limit as i? — )■ oo, see Corollarv 13.11 /^^(o-J?) W '^^ 
a limit, that we denote by /i. 

Let y = kx. On one hand, we have \y\^~^V'ip{y) = \x\^~^'^'il)^{x) and on the other, 
\y\^~^\/ip*[y) = \x\^ ip* [x) . Therefore, taking k = kn and \x\ = 1, which implies |y| = k^, 
and using (j3.7p . we get 



^IVipiy) - Vip*{y) \ = \Vip''"{x) - Vip*{x)\ <e if |y| = kn and n is large. 



Hence, 



/ ^dS-C*{2-N)n''-'(dBiO;l)) = [ 



dip dip* 
dr 



dS 



< Ce 



if n is large. We conclude that, C* = fi ((2 - N) n^~^{dB{0; 1))) \ 

Rephrasing the limit. Take y = kx and |a;| = 1. In the above steps we have proved that 

\y\^-'^+'Df^iP{y) =k^~'^+'D>^iP{kx) = D>^ip''{x) ^ D^iP*{x) as |y| = A: ^ oo. 

Since D^ip*{x) = \y\^^'^'^^^^D^ip* {y) (remember that \x\ = 1), this is just ()3.6p . □ 
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Remark 3.3. To get (|3.6p it is enough to have J G (7max(3,s)_ -^yg need at least J G C^, no matter 
the order of the derivative we are considering, in order for A-^ to be integrable at infinity, which 
is one of the main ingredients of the proof. 



4. Conservation law and asymptotic mass 

We start our study of the evolutionary problem (jl.ip with a conservation law that will be of 
the most importance in order to establish the behavior of the solution in the far field. 

Proposition 4.1. Let u be the solution to (jl.ip and <j) the solution to (jl.Sp . Then, 



(4.1) 



u{x, t)(j){x) dx 



Uo(x)(/)(x) dx := M* for every t > 0. 



Proof. Since u{x, t) = 4>{x) = for x £ T-L, we have 

d 
di 



u{x, t)(j){x) dx 



I / 



ut{x, t)(p{x) dx 

J{x — y) [u{y, t) — u{x, t))(/)(x) dy dx 
J{x — y) [4>{y) — 0(x)) u{x, t) dy dx 



□ 

Though mass is not conserved, there is a non-trivial asymptotic mass, which coincides with 
the conserved quantity M* . 

Proposition 4.2. Let u he the solution to (jl.ip and M* = J^pf uq{x)(I){x) dx. The mass of the 
solution at time t, M{t) = / u{x,t)dx, satisfies M{t) — )• M* as t ^ oo. 



Proof. In fact. 



u(x, t) dx - M* 



u{x, t) (l — ipix)) dx 



< I u{x,t)dx+ / u{x,t)(^l — cf){x)^ dx 

'\x\<riVi J\x\>riVt. 

< Cr?^ + C7||no||Li ^ 



<Ct-^/^{r]Vi)'' + C- 



ivVt) 



N-2 



Hence, choosing rj so that Crj^ < e and taking limsup^^g^ we get. 



lim sup 



u{x,t) dx - M* 



< £. 



We conclude by letting e — 0. 



□ 
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5. Far-field limit 

In this section we study the large time behavior of the solution to (jl.ip in sets of the form 
> 5t}, with 6>0. 

Theorem 5.1. Let u be the solution of (jl.ip and M* the asymptotic mass given in ()4.1|) . Then, 
for every 5 > 0, 

(5.1) t^/^||u(x,t) - M*ra(x,t)||ioo({|^|2>5f}) as t ^ oo. 

Proof. Let to > and t > 2tQ. From the representation formula (|2.9p we have 

t^^^\u{x,t) - M*Ta{x,t)\ < t^/h^^'~''''^u{x,to) 



.N/2 



ui{x -y,t- to)u{y,to)dy - M*ra{x,t) 



l2 



+ r e-(*-^);f-H(x)(J * u{-, s)){x) ds 

Jtn 



h 



.N/2 



to JR'^ 



uj{x -y,t- s)Xu{y){J * s)){y) dy ds . 



h 

Let 6 > and |xp > 6t. We have the following estimates for the quantities If. 
— )■ as t — )• oo. 

► The estimate for I2 is the most involved, and is split in three parts, 

h < t^^^ / 1^(2; -y,t-to)- Ta{x -y,t- to)\ u{y,to) dy 



h.i 



Taix-y,t- to)u{y, to) dy - M{to)Taix, t) 



h,2 



+ t^/'^r^{x,t)\M{to)-M* 



Using (pTUjl . 

h,! < t^/'^\\u{-,to)\\Li(j^N)\\u{-,t - to) - Ta{-,t- to)\\Lo 
= t'^/'ho||Li(RiV)o((t-to)"^). 



Thus, limt^oo h,! = 0. 

The well-known asymptotics for (local) caloric functions imply limt-5.00 -^2,2 = 0. 
By Proposition K2\ if to is large enough, 12,3 < Ca\M{to) - M*\ < e. 
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► For to large enough, x ^ Ti. Hence, I3 = 0. 

► If tQ is large enough and y & Ti, \x — y\ > \x\/2. On the other hand, by ()2.5p < u{x,s) < 
Cs-^/^ and therefore, < (J * «(•, s))(x) < Cs-^/^. Thus, using the space-time bound ()2.14p 
for w, 

i> < ct-'^ f / < c(^) < crf-<--« < . 

J J 'hi \ y \ II 
if to is large enough. 

Gathering all the above estimates, we get 

limsup \\t''/\{x,t) - t^/2jvrr„(x,t)||^^(|^p>5^) < 2e. 

We finally let e ^ 0. □ 

As a corollary of the outer behavior, we can do better in the proof of Proposition 14. 21 to obtain 
the rate of decay of the mass to its asymptotic limit. 



Corollary 5.1. Under the assumptions of Proposition we have 

JV-2 JV-2 

(5.2) M{t) = M* + Kt — — + o{t — —), 

where K = C*M* J^r, Ua{0\C\'^~^ d^, C* = lim (1 - (t>{x))\x\^-'^ . 

Proof. We have from the conservation law ()4.ip that 

JV-2 iV— 2 f 

t—{M{t)-M*)=t— / u{x,t){l- <l){x))dx. 
In order to estimate this last integral we split it into two parts, 

JV-2 f 

t 2 / u{x,t){l - 4>{x)) dx = 

JV— 2 f JV— 2 f 

t 2 / u{x,t){l - (/){x)) dx + t 2 / u{x,t){l - (l){x)) dx 

J{|a;|<(5il/2} J{|a:|><5ti/2} 

^ V ' V ' 

h h 

ioT 6 > and t big enough. Using the estimate < 1 — <l){x) < ci|xp^^ and the bound 
< u{x, t) < Ct~^/'^ we get < /i < 05"^ . For we perform the change of variables x = ^t^^^, 
and get 



= C*M* [ Uaimi^^^'dC. 
J||fl><5l 



Hence 



C*M* [ C/„(C)|CP~^c^? < liminftV(M(t) -M*) 

< limsupt^(M(t) - M*) < C5^ + CM* [ UaiOl^l"^'^ d^, 



'm>s} 

from where the result follows just by letting (5—7-0. □ 
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Remark 5.1. The amount of mass, Ml{u), lost in the evolution is given by 

Ml{u) := / uo{x) dx — lim / u{x,t)dx= / {1 — (j){x))uQ{x) dx > 0. 

Therefore, the influence of the hole structure is felt at the asymptotic level through the projection 
of the initial data onto ^ := 1 — <j), which represents in this way the dissipation capacity of Ti. 

Notice that ip is the L-harmonic function defined in Q that takes value 1 on ?^ and at 
infinity. Hence, by analogy with the standard (local) potential theory, we may say that is the 
function that measures the L-capacity of Ti, by means of the formula 

cap(«)= mt n Ji.-y)(ui.)-u(y)fi.iy. 

{u>i onnyJn Jn 



6. Near-field limit 

In view of Theorem 15. H what is left to complete the proof of Theorem 11.11 is to show that the 
limit ()1.6p is valid uniformly in sets of the form < 6t} for some 6 > 0. This is the goal of 
this section. 

Theorem 6.1. Let u be the solution to (jl.ip . There exists a value 5 > such that 

(6.1) \imt''/^uix,t) - MXx)r„(x,t)||ioo(||,|2<5j|) = 0. 

Instead of ()6.ip . we will prove 

(6.2) lim t^/^\\u{x,t) - M*0(x)w(x,t)||ioo(||,|2<5t|) = 0, 

which turns out to be equivalent, thanks to (j2.10p with |/3| = 0. Notice also that, using 
again ()2.10p . together with the asymptotic behavior of (p, we may rewrite ()5.ip in a similar 
way, 

(6.3) lim t^/^\\u{x,t) - M*(l){x)uj{x,t)\\Lo.^{\^\2>st}) = 0. 

In order to prove ()6.2p we will construct suitable barriers approaching the asymptotic limit 
as t goes to infinity. We choose k G (0, min{l, — 2}), 7 G (0, (A^ — 2 — k)/2), and then define, 
for any K± > 0, 

JV + K 

(6.4) v±{x, t) = (p{x)uj{x, t) ± K±t ~z{x), 

with z as in Lemma 13.21 z(x) = {\x\'^ + b)~^ . We take b large enough, so that (|3.ip is satisfied. 

Lemma 6.1. Let R > and v± as above. There exists a value (5* > 0, independent of K±, such 
that for all 5 G (0, 6*), K± > 1, we have 

dtv+ - Lv+ > 0, dtv- - Lv^ < 0, R^ < \x\^ < 6t. 



Proof. On one hand. 



dtv+{x,t) = (j){x)dtuj{x,t) 



N + K. _ JV+K 

K+t 2 z{x) 

2t ^ ^ ^ 



NONLOCAL DIFFUSION IN DOMAINS WITH HOLES 19 
On the other hand, using that (p is L-harmonic and (|3.ip . we get 

Lv+{x,t) <(l){x)Luj{x,t) + / J{x - y){(l){y) - (j){x)){u}{y,t) - uj{x,t)) dy 

Hence, since to solves ()2.13p and > 0, 

dtv+ - Lv+ > - / J{x - y)|(A(y) - (l){x)\\u}{y,t) -uj{x,t)\ dy 



A 



N+. ^ /2a-y(N -2-2-f) N + k 
+ K+t — —z{x)' " 



|x|2 + 6 2t 

V 



Estimate for A. We first notice that there is a constant C > such that 

\(k{x) - 4>{y)\ < if|x-y|<l. 

Indeed, since \4){x) — 4){y)\ < 1, this is obviously true for |x| < 2, with C = 2^~^ . On the other 
hand, using Taylor's expansion, 

(6.5) |</)(x) — < max max \D^(j){^)\ ii\x — y\<\. 

l/9|=U6B{x;l) 

The result for |x| > 2 now follows from the estimate (j3.4p with |/3| = 1, since |^| > \x\/2 for all 
^ £ B{x] 1) whenever |x| > 2. 

Combining (|6.5p with (|2.1ip . we conclude that there is a constant C > such that 
A < , , = , ' < — —TT-r if xP < 5t. 



Estimate for B. If < \x\'^ < 6t and 6 < "[^^^"^g^^g' , we 

N + K S{N + k) _ 6{N + k){R^ + b) 6{N + K){R^ + b) aj{N - 2 - 2^) 
2t - ~ R'^\x\'^ + blxl"^ - i22(|x|2 + 6) - |x|2 + 6 ' 

Thus, with this choice of 6, 

^ ^ a^jN - 2 - 27) 
The above estimates for A and B yield, if > 1, 

On the other hand, since 2(7 + 1) < — k, 

(|x|2 + b)T+i (i;2 + 5)7+1 

Therefore, taking 5 small enough, we get 

dtv+ - Lv+ > 0, R^ < \x\^ < 6t 



have 
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for all Kj^ > 1. Notice that the threshold value of 6 does not depend on Kj^. 

An analogous argument leads to the statement concerning since e"*J(x) < Ct'~^\x\^~~^ . 

□ 

Proof of Theorem \6.1\ Let -R > such that -6(0; R) C Ti. As a consequence of Lemma l6.ll we 
know that there is a value 5 > such that 

dtv+ - Lv+ > 0, X E O n < 6t}, t > 0. 

Let now e > 0. By (|6.3|) . (|'2.10p with |/3| = 0, and the asymptotic behavior of (j), there exists 
to > such that u{x,t) < {1 + e)M*(j){x)uj{x,t) if ((Jt)^/^ < [^.j < (5t)i/2 + l, t > to. Hence, for 
any > 0, 

u{x,t) <{l + e)M*v+{x,t) if {5tf/^ < \x\ < {5t)^/^ + 1, t > to- 

On the other hand, it is obvious that there exists > 1 such that to) < {l + e)M*v^{x, to) 
if |x| < ((5to)^/2 + 1. Finally, u{x,t) = < (1 + e)M*v+{x,t), x e U, t > t^. Therefore, by 
comparison, see Theorem 12.31 and Remark 12. 3| 

u{x, t) < (1 + e)M*v+{x, t), \x\ < ((5t)^/^ t > to. 

Hence, using the decay estimate ()2.1ip . if |x| < ((5t)-^/^, t > to, we have 

t^/^ {u{x,t) - M*<)){x)io{x,t)) < eM*(/){x)t^/^uj{x,t) + {l + e)M*K+t-^z{x) 

< eM*C + {l + e)M*K+t-^b~'-'. 

Letting t — ?• cxo and then e — t- 0, we conclude that limsupt^^^ (u(x,t) — M*(j){x)Lo{x,t)) < 0. 

An analogous argument shows that 

(1 - e)M*v_{x,t) < u{x,t), \x\ < ((5t)^/^ t > tg. 
Hence, liminf t^/^ [u{x,t) - M*(j){x)io{x,t)) > 0. □ 

t—^co 

As a corollary, we obtain the behavior on compact sets. 

Corollary 6.1. Let u be the solution to (jl.ip . Then, 

Ad* 

(6.6) t^^'^u{x,t) — )• — uniformly in compact subsets of] 



Remark 6.1. The limit ()6.6p holds uniformly in sets of the form Q, n {|xp < 6(t)t}, as long as 
lim 6{t) = 0. 
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